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ABSTRACT 


The  determination  of  an  orbit  by  using  two  heciocentric 
positions  of  a  celestial  body  is  a  boundary-value  problem  of 
mathematical  physics,  and  may  be  reduced  to  the  solution  of 
a  nonlinear  integral  equation.  Methods  of  successive  approx¬ 
imation  to  the  solution  are  discussed,  and  convergence 
conditions  and  error  estimates  are  derived. 

The  problems  of  perturbed  motion  and  the  determination 
of  orbits  based  on  three  geocentric  observations  are  also 


considered. 


ON  THE  DETERMINATION  OF  ORBITS  BY  THE  SOLUTION 


OF  A  SYSTEM  OF  INTEGRAL  EQUATIONS 
V.  T.  Gontkovskaja 

Bjull.  Inst.  Theoret.  Astonom.  £(  1962),  283  -  29  8 
Translated  from  Russian  by  L.  B.  Rail 

One  of  the  basic  problems  of  celestial  machanics  -  the  determination  of 
an  orbit  from  two  heliocentric  positions  -  may  be  considered  to  be  a  boundary- 
value  problem  of  mathematical  physics.  Until  recently,  such  problems  have 
been  solved  by  integrating  differential  equations  with  given  initial  conditions. 
This  reduction  *rom  one  problem  to  the  other  turns  out  to  be  possible  only 
because  the  perturbed  motion  differs  only  slightly  from  Keplerian  motion  in  the 
two-body  problem,  the  general  solution  of  which  is  known  for  any  statement  of 
it  (either  as  an  initial-value  or  a  boundary-value  problem). 

At  the  present  time,  a  whole  series  of  problems  has  arisen  in  connection 
with  the  conquest  of  space  in  which  the  perturbations  are  very  large,  and  the 
Keperlian  orbit  cannot  serve  as  an  approximate  solution.  Thus,  the  determination 
of  an  orbit  in  such  cases  is  only  possible  directly  from  the  boundary-value 
problem,  for  example,  by  reducing  it  to  an  integral  equation. 

This  method  was  first  proposed  by  Bucerius  in  his  fundamental  paper, 
"Bahnbestimmung  als  Fla ndwert problem.  "  Bucerius  not  only  stated  the  problem, 

but  also  proposed  a  method  for  its  solution.  However,  the  method  of  successive 
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approximation  proposed  by  him  for  the  solution  of  the  system  of  integral  equations, 
to  which  the  problem  of  determination  of  orbits  from  two  heliocentric  positions 
reduces,  contains  a  whole  series  of  inaccuracies  and  does  not  have  a  rigorous 
foundation.  Indeed,  for  the  method  of  successive  approximation  proposed  by 
Bucerius  for  the  solution  of  systems  of  integral  equations  differing  somewhat 
from  those  given  here,  some  conditions  given  for  convergence  are  obtained  as 
a  result  of  non-rigorous  reasoning,  and,  what  is  very  essential,  Bucerius  does  not 
have  any  error  estimates  whatever  for  his  method. 

In  this  paper,  two  methods  of  successive  approximation  are  considered  which 
permit  the  solution  of  the  problem  posed  by  Bucerius.  Sufficient  conditions  are 
obtained  for  the  convergence  of  the  sequence  of  approximations  to  the  solution, 
the  uniqueness  of  the  solution  obtained  is  proved,  and,  in  addition,  estimates 
are  given  for  the  error  committed  by  terminating  the  calculation  at  the  nth 
approximation. 

§1.  Statement  of  the  Problem 

The  differential  equations  of  motion  in  vector  form  are: 

2 

F(t)  =  -  k  ( 1  +  m) - —  .  (1) 

Irltll3 

We  introduce  the  vector  r^f  t) ,  which  satisfies  the  equation  of  uniform  rectilinear 
motion 

U  l>  ■  0  • 

in  the  form 
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where  and  are  determined  by  two  heliocentric  positions  at  the  times  tj 
and  t, .  Thus,  the  vector  r(t)  -  r„(  t)  satisfies  a  second-order  differential 
equation  with  zero  boundary  conditions,  which  properly  completely  define  a 


solution. 


By  means  of  the  Green's  function 


G<t,t')  = 


(t-tjM^-t*) 

h-h 


t<t* , 


G{  t,  t')  =  G(t*,t) 


the  problem  may  be  reduced  to  the  solution  of  the  nonlinear  integral  equation 


r(t)  =r0(t)  +  k2(l+m)  /  G(t,  t') 


r(t') 

r(t»)i; 


If  the  time  t^  is  taken  to  be  the  beginning  of  the  measurement  of  time, 
and  the  difference  t^  -  t^  as  a  unit  of  time,  then  the  time  t  varies  over  the 
interval  [0,1],  The  kernel  of  the  integral  equation  takes  the  form 


K(t,t«)  = 


t(l-t'),  t<t' 


t'(l-t),  t>f 


and  equation  (  2)  may  be  written  in  the  following  fashion 


r(t') 


£(t)=r  (t)  +  T  J  K(  t,  t' )  — - -j-dt'  , 

0  lr(  t')l 
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where 

t  =  k^l  +  m  ( t2  -  tj)  . 

Since  the  vectors  r  f X  ^  »X  2  are  c°P^anar»  relationship 

r(t)  =  n](t)r1  +  n2(t)  r2  ,  (4) 


holds,  where 

n1(t)  =  n1(t,t1,t2,r1,r2)  , 
n2(t)  "n2(t,  tj,  t2>£  1  »X  2)  * 

Interpreted  geometrically,  npn2  are  the  ratios  of  the  areas  of  the  triangles 


n.(  t)  = 


n,(  t)  = 


Uill 


l'fc/  [r.  r  ]  »  2'  '  [  r  r  ]  * 

~  1~  2  ~  1~  2 


With  the  help  of  relationship  (4),  the  integral  equation  (4)  reduced  to  the  two 
integral  equations 

1 


n(t)  sl-t-l-r  /  K(t,t') - - - - 

0  |n1(t')r1+n2(t*)r2l 


dt' 


n  (t)  »t  +  r 2  /  K(t,t')- 


n2(t') 


n1(t')r1+n2(t,)r2 


dt' 


>  (5) 


ni<  t' ) £  i+  n2<  t' ) L  2 1 ' 3  =  I  nf(  t' )  rf  +  n2(  t* )  r2  +  2ni(  t' )  n2<  t* )  ( £  1 X  2)  1 3  ^ 1 1 


§2.  A  Method  of  Successive  Approximations  Applied  to  the  Nonlinear 
Integral  Equation 

To  solve  equation  (  3) ,  Bucerius  (  1950)  proposed  an  interative  method  which 
consisted  of  a  sequence  of  applications  of  the  mean-value  theorem  to  the  integral 
on  the  right  side. 
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Suppose  that  rQ(  t)  is  taken  as  the  initial  approximation  to  the  solution  of 
the  integral  equation  (  and  the  initial  approximations  to  the  ratios  of  the  area  as 


T1  7  2 

—  ,  — ) ,  then  the  first  approximation  is 


1 


2  1 


Ij<  t)  =rQ(  t)  +  r^/K(  t,  t')-^°-  -jdtf  -rQ(  t )  +  ~  / K(t,f)  rQ(f)  df 


0  lr0(f)l 


rl  0 


The  average  value  of  the  quantity  |r  ( t')|  in  the  interval  0  <  t*  <  1  is 


denoted  by  . 


Introducing  the  notation 


Vt  - 

ri 


then 


£,( ')  =£0(  t)  +  xjK(  t,  t')r0(  t')df  . 


The  second  approximation  has  the  form 


2  } 

r  <  t)  =r  't*  +  t  / K(t,f)-  df  . 

0  |r,(t')|J 


If  the  expression  for  r .( t')  is  substituted  into  it,  then 

1  1 

r  2<  t)  =r0<  t)  +  x2  / K<  t,  f)  r0(  f)  dt'+X^  /  K2(  t,  f>rQ(  f)  df  , 


where 


> 


and  Kg(  t,  t1)  is  the  second  iterated  kernel. 
For  the  vth  iteration,  one  will  obtain 
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1  1 

Xv<  t>  -I  <t)  +  \,/K<t,f)r  (t'lcU'  +  X  X  ,/K  (',*')  r  ( •'  >  “f  +  ...  + 

0  0 
1 

+  Vv-l---Xl/Kv(t>t,,io<t',dt'  >  (6) 

0 

2 

x  =  "T  » 

v  3  ’ 

r 

v 


K  ( t,  t* )  is  the  vth  iterated  kernel. 
v 

2 

All  \  satisfy  the  inequality  \  <\=~,  where  h  =  min  r(  t)  (  0  <  t  <  1) 

V  y  — “  J  ~  —  — 

h 

is  the  perpendicular  issuing  from  the  origin  of  coordinates  to  the  chord  r .  r_  . 

2  ~1~2 

Letting  v  go  to  infinity  and  replacing  all  X  in  (  6)  by  — r  ,  one  may  consider 

V  h 

the  series  obtained  to  be  the  Neumann  series  for  the  equation 

2  1 

r(t)=r0(t)  +  ^y  /  K(  t,  t')  r(  t'Jdt*  . 


A  lower  bound  for  the  radius  of  convergence  of  this  series  is  (Wiarda,  1930) 

I  1  1  2  \-1/2  __ 

X  <  /  /  [  K(  t,  t*)j  dtdt'  =  3-s/ 10  , 

\  0  0 


or  x  <X  ,  where  Xj  is  the  first  characteristic  value  of  the  kernel  K(  t,  t'). 
Consequently,  if  all  X^  <  X^ ,  then  the  series  (  6)  converges. 

In  applying  the  iterative  process  described  above  to  the  solution  of  equation 
( 3),  Bucerius,  generally  speaking,  started  from  the  mistaken  assumption  that  the 
X^  do  not  depend  on  time,  in  the  sense  that  they  are  the  same  for  each  orbit. 
Indeed,  X^  =  X^(  t),  and  one  has  to  prove  the  assertion  that  all  ry  ^h,  from 
which  the  convergence  of  the  process  of  successive  approximation  follows. 
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§3.  Solution  of  the  Nonlinear  Integral  Equation  by  the  Method  of 
Successive  Approximations. 

We  apply  to  equation  ( 3)  the  method  of  successive  approximations  in  the 
following  form. 

We  first  take  r(t)  =£0(t)  and  define  r^t)  by  the  relationship 

2  rl  l0iV) 

~l(t)  =~0(t)  +T  /K(t»t')  3 - dt’  • 

0  rj(f) 

Continuing  in  the  same  fashion,  we  obtain  an  infinite  sequence  of  functions 

J0(  *)>•••  >£n(  *)»•••  > 


which  satisfy  tl  e  recurrence  relations 


~rn(t)  =~ro(t)  +  T  /«*»**> 

0 


1  r  .(f) 

... .  ...~n-l 


df  . 


(7) 


00  converges  to  a 


We  shall  show  that  the  sequence  of  vectors  lr  ( t)  } 

~n  n=0 

function  which  is  a  solution  of  equation  ( 3). 

For  the  proof,  consider  the  difference  of  two  successive  approximations 


•in' 11  'Vl*  ""Vwt.  »*» 
0 


-rn-l(tl)  -rn-2(t'> 


dt'  . 


The  difference  enclosed  in  square  brackets  may  be  written  in  the  following 


manner. 


r3  ,(t')  r3 
.  n-1  n-1  -i 


*n-2(t'>  -rn-2(tl) 

v.(t,»’v2<t'b 


r  .( t')  -  r  ,( t') 
~  n-1  ~  n-2 


r  ,(t')-r  (f) 

n-2  n-i 


+  r  ( f )  — - ~ ■ — 

+  ~n-2' 1  r  .(f)  r  (f) 
n-i  n-2 


1 


1 


1 
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In  so  far  as 


J  K(  t,  t1)  dt*  <tj>  for  0<t<l, 


'~n-l~~n-2 


n-1 


1 


r  rn-2~rn-l  1  I  _ 

~n-2  r  ,r  ,  '  2  r  ,r  _ +  2 

n-1  n-2  r  ,  n-1  n-2  r  .. 

n-1  n-2 


and  since  the  difference  of  two  sides  of  a  triangle  is  less  than  the  third  and 
r 

|— — - 1  £1,  one  has  that 

n-2 


,  max  |  r  ,(t)-r  ,(t)| 

2  «  ^  .  ~  n-1  ~  n-2 

,3  3. 

0  < t  <1  n-rrn-2* 


(8) 


or 


( 


0<t<l 


n-1 


2  min  r, 
0<t  <  1  * 


(9) 


k  =  0  ,  1 ,  2 ,  . . .  ,  n-1  . 


Setting 


r  ( t)  -  r  ( t)  =  A  ( t) 
~  n  ~  n-1  ~  n-1 


> 


we  note  that 


n-1 


£n< [)  "£0( ')•»[£,(« -£0<  «!  +  [j:2(t)-£l(  «]+••• +  t£n,t)-rn_i(  91  =r0(t)  +  ^  Afc(  t) 

k=0 


and 


r„<  0<ro<  t)+A0(  t| 


2 

1  +  -T  + 

1 

(  .  ' 
T 

n-1 

1  +  3  + 

2  min  i  (  t) 
n 

L  0  <t  <i 

3 

2  min  r  ( t) 

Vo<t<l  n  J 
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since  ( 10)  is  valid  for  all  r  ( t),  the  limit  of  which  is  R(  t) . 

~  n  ~ 

The  function  R(t)  satisfies  equation  (3).  In  fact,  by  virtue  of  the 
uniform  convergence  of  the  sequence  {r^f  t) }  to  R(  t),  by  letting  n  -*  «  in 
formula  (  7) ,  one  obtains  that 

2  rl  5(t'> 

R(t>  =r  ( t)  +  T  - dt'  . 

0  R(t') 


We  shall  prove  that  R(t)  is  the  unique  solution  of  equation  (3)  in  the 

region 

2 

- -  2  <1  •  (12) 

2  min  r  ( t) 

0<t  <1 
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Indeed,  suppose  that  there  exists  another  solution  £(  t)  which  satisfies 
condition  ( 12) .  Then,  by  performing  a  calculation  similar  to  the  preceding,  one 
obtains  that 

Bit'*  p<  t»)  4|p(t')-R(t')| 

R3(t')  p3(  t')  min  [  R3(  t1) ,  p3(  t')J 

0<t'  <1 


Consequently, 


!  R(  t)  -C(  t)| 


t2/k(  t,  t') 
0 


g(t') 
R?(  f) 


p(t') 
p3(  t') 


t2  I  R(  t)  -£  ( t)  | 

2  min  [  R3(  t)  ,P3(t)] 
0<  t<  1 


(0<t<l), 


or 


I  r(  t)  -£(  t)i 


/ 


min  l  R  ( t) ,  P3(  t)] 
0<t<l 


<  0  . 


The  quantity  in  parentheses  is  positive,  by  virtue  of  the  satisfaction  of 
conditions  ( 11)  and  ( 12) .  Dividing  both  sides  of  the  above  inequality  by  it,  one 
obtains  that 

I R ( t)  —  p { t) |  <o,  ( o  <  t  <  1)  , 

from  which  it  follows  that 

I  R  ( t)  -  p  ( t)  |  0  ,  (0  <  t  <  1)  . 


We  now  derive  a  condition  which  guarantees  the  satisfaction  of  inequality 

(10). 

It  follows  from  equation  (  7)  that 


r„<  >  r0<  t> 


2 


8  min  r2  ( t)  ] 
0  <t<  1 


(13) 
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It  is  easy  to  seen  that  rQ(  t)  is  larger  than  the  altitude  h  from  the  origin 
of  coordinates  to  the  chord  joining  the  points  ( x^,  y^,  z^)  and  (  x2,  y 2,  z2) . 


Then, 


min 

°<t<l 


r  ( t)  >  h  - 
n  — 


2 

T 


8  min  r 
0<t<l 


2 

n-1 


(t) 


(14) 


Assume  that 


min  r  (t)>kh,  (  m  =  0, 1, . . . ,  n  - 1)  ,  (15) 

o<t  <1  m 


where  k  <  1  is  a  positive  number  which  we  shall  define  below.  Note  that  this 

inequality  is  valid  for  m  =  0  .  We  assume  that  it  is  valid  for  all  m  to  n-1, 

inclusive.  With  the  help  of  inequalities  ( 14)  and  ( 15) ,  we  obtain  that 

2 

min  r  ( t)  >  h  - — y-r-  . 

0<t<l  8k  h 


Requiring  that 


2 

h  --T-r  >kh  , 

8k  h 


(16) 


we  obtain  that  min  r  ( t)  kh,  so  that  we  have  proved  this  inequality  for  all 
°<t<l  n 

n  by  means  of  the  principle  of  mathematical  induction.  Thus,  by  virtue  of  ( 10) , 
the  inequality 


2k3h 


3  <  1  } 


(17) 


is  a  sufficient  condition  for  convergence. 

4 

Both  inequalities  ( 16)  and  ( 17)  are  satisfied  if  one  chooses  k  <  —  . 

—  5 

Actually,  ( 16)  is  a  consequence  of  inequality  ( 17)  if  one  requires  further 


that 
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from  which 


h  - 


2kV 

8k2h2 


>kh  , 


1  -Tk >  k 
4  — 


> 


or 


k 


Consequently,  the  process  of  successive  approximations  described  above 
converges  to  a  solution  of  our  integral  equation  in  case  that 


.3  .  125  2 

h  -128  T  ' 


for  which 


min  r  ( t)  >  i  h  . 
0<t<l  n  5 


We  now  obtain  an  estimate  for  the  error  committed  if  the  calculation  is 
halted  at  the  nth  approximation: 


where 


Since 


00 


k=n 


n 

— - a 

1-  a  O’ 


a 


2 

T 


2  min  r 
0<  t<  1 


3 

n 


(t) 


<  !£L  I 

-  128  . 
n 


2 

T  • 


/  K<  t,  f ) 


Jo<  t~) 
ro(t'> 


dt' 


2 

_ T _ 

8  min  r2(  t) 
0<t<l 


9 


finally, 
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\bui 


|r(t)-r„(t)|  <j^-h 


/  2  1 

n 

125  t 

128  ,3 

'  h  ; 

125  J_ 

1_ 128  ,3 
n 


§4.  The  Application  of  Newton's  Method  to  the  Solution  of  a  System  of 
Nonlinear  Integral  Equations. 

Newton's  method,  sometimes  called  the  method  of  tangents,  is  one  of  the 
most  effective  methods  for  obtaining  solutions  of  algebraic  equations  in  the  case 
that  an  initial  approximation  is  known  for  the  value  of  the  solution.  The  successive 
approximations  are  defined  by  a  formula  of  the  form 


x  .  =  x  - 


f(xn) 


n+1  n  f*(  x  ) 
n 

This  method  extends  to  systems  of  algebraic  equations.  It  may  be  applied 
(  KantoroviS,  1949)  to  classes  of  arbitrary  nonlinear  equations,  in  particular,  to 
classes  of  nonlinear  integral  equations. 

We  shall  apply  Newton's  method  to  the  system  of  nonlinear  integral  equations 

1 


fc)  -1-72/  t,  t',  t') ,  n2(  t'))dt'  , 

n^t)  =  t) -i-T2/  U2(t,  t',  n^  t') ,  n2(  t'))dt'  , 


(18) 


for  which  initial  approximations  to  the  solutions  may  be  taken  to  be  n^ ( t) ,  n2(  t) . 
To  obtain  the  first  approximations,  one  adds  to  n?(t)  ,  n?(  t)  the 


corrections 


A2  =  llj ( t)  -  Jljc  0  , 

a2  =  4( i|  -  4(  i>  - 
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determined  from  the  linear  integral  equations 


?  a}+ 

8Vo  1 


air)  *2}^'  =  j2/  ui(  t,  t*f  *1  ®c  t") ,  n°(  t*))  dt«  , 
2  0  ^  0 


!-tY( 

i-T7  {(^]0AV(^)oa2}dt'  =  T7U2(t’t,’-l<t,,’-2<t,,,dt  ’ 


(19) 


2  2 

The  second  approximations  n^(  t) ,  n  t)  will  be  obtained  by  determining 
2  2 

A^  ,  A^  from  the  system 


?  7  IT/ 0Ui  ?  /9U,  \  n  1  7  1  1  1 

1-T  /\.l'9^l1Al+l'Sr2)1A2jdt'  -l“)'ILl(t)+T  /U1<t,t',nl(t'),n2(t'))df  , 


!  ZflfldUZ 

•■-r  /  Ik 


* 


au. 

_ t 

an,  I 


2  '1 


^jdt'  =  n°(  t)  -n*(  t)  +T2Ju2(t,  t',  nj(  t'),n*(  t'))dt' 


with 


n  *(  t) 
*2(t) 


=  n[(  t)  +  A2  , 
=  I>2<t)+  A2  , 


finally,  the  nth  approximation  is  determined  by  solving  the  following  system  of 
linear  integral  equations: 


An  2 

VT 


0 


au 


1 


a"  + 


^Vn-l  1 


au 


1 


.  n 


n-1 


=  n  ”( t)  -  n  "J(  t)  +  t‘J  U1(t,t>,n"',(t'),  n2“(f))dt-  , 


A"  )df  = 


n-1. 


n-1. 


.  n  2 

VT 


au. 


v  3n1 


n-1 


,n 

V 


au. 


>  (20) 


n 


1  9n2  'n-1  lJ 


*  n“(  t) -n"_1(  t>  +  T2/u2(t,t',n"‘1(t>),  njj'Wdt' 


dt*  = 
n-1. 


J 
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**  » rm,  «r-  ■-»  j».« 


In  the  case  of  interest  to  us,  system  (  5),  the  functions  U  and  the 

9U.  au2  1 

derivatives  - —  ,  — —  { i  =  1,  2)  have  the  following  form 
9n,  ’  9n,  ’  ' 


U,  = 


K(  t,  t*)  n](  t«) 


1  r  2  2~  2  .  .,3/2  » 

[  n!rl  +  n2  +  2nin2(~l~2)] 

K<  t,  t')  n  ( t») 


U2  ~r  2  2  2  2 


[  nl  r  1+  n2r2  +  Zn^n^  r  1  r  2 ) ] 


3/2 


9UX  3K(  t,  t')  nj  t')  [  n2r2  +  nj(  r  1  r  2)] 
9n„  5 


aui  k<  t,  t')  r,  3nitniri+n2<~ri-rz>n 

1  r3  v-  r2  J 


9n 


9U2  3K(  t,  t*)  n2(  t*)  [  n^j  ] 

7  ’ 

K(t.f,r.  3n2ln2r2tnl(-r|^2|l  1 

3n2  r3  l  r2  J  *  J 


(21) 


moreover. 


max  (n  ( t) ,  n  ( t)  }  ''j 

|/u.(t,t')dt'|±j^ - 3 -  , 

0  min  r  ( t) 

CK  t<l 

la  M2  max  (n2(  t) ,  n2(  t) } 

i  i  „  1  0<t<l  \ 

— dt'll-  - Z -  ,  > 


0  3ni 


-  8 


min  r  ( t) 

°<t<_l 

M2  max  (n2(  t) ,  n2(  t) } 


0  3ni 


min  r  ( t) 
0<  t<  1 


J 


(22) 
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vt'  V*  >**/*•*  • 


2 

where  M  denotes  the  greatest  of  the  four  following  quantities: 

2rl+r2  +  I(^i£2)I  » 

rl+  ll£1i2)l  , 
rl  +  l(£,£2)l, 
r2+2r2  +  l<il ;-r2)l- 

We  shall  obtain  the  conditions  for  convergence  of  the  process  of  successive 
approximations  defined  by  formulas  (  20).  We  write  them  In  the  following  form: 

n"(t)  =  n°(t>  +  t2  /  Uj(  t,  f,  n"'1(t'),  n"*'(f))df  +  1 


,  1(  9U  fdu  l  *s 

■/{■ST  At1  An)df 

0  U8"l  /„_!  1  1  8n2in.!  2jdt 


n2<  t>  ="2(«  +  t2/  U2(t,f,nn-1(f),n"-1(f))df  + 

0  1  c 

7  Ifld  u_  fdu 

+  T  n  Ur  Ain+  ir  a"1,.  . 
0  dnl  ;n-l  1  an2  'n-l  2J 


(23) 


Subtracung  n^(  t,  from  nj't  t,  and  n^(  t,  from  njt),  one  obtains 
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>  N  N 


a"  =  T2/ [  Uj(  t,  f,  njl_1(  t'),  n^"1*  t‘>)  -  Uj(  t,  f,  n"'2(  t>),  nn2~2{  f))]df  - 


-r2/rp-  a"-‘+(9A]  A"-Adt.  +  T2/&)  A?+f  — ^  I  Aj'ldf, 

ol'8nl'n-2  i  8n2  ln-2  2  J  oli8nl'n-l  1  ^8  2'n-l  2J 


Aj  =  t2/  [  U2<  t,  t\  n"_1(  f),  n2‘‘<  f))  -U2(  t,  t1,  n"'£(  f),  f»]dt*  - 


n-L  lv  n-l. 


2  fVfaU2  \  An.l  l,u2|  n-ll  ,  2,y/8U2|  n  J8U2  )  .n\ 

-7  f  it  Ai  +  sd  A  jde+Ti\*r,  ,Ai  +  \ir„  ,a2  r 

0  k'  1  >n-2  *  2'n-2  J  0  V  l'n-1  v  2'n-l  J 


Using  the  notation 


An  =  max  {)  aJ1!  ,  |  A^|  } , 
0<  t  <  1 


n  /  n  n> 

n  =  max  in.  ,  n  }  , 

0<t<  1  1 


and  representing  the  difference 

n2_1(  t>»  -  U,<  t,  t>,  n"’2(f),  n"_2<  t>)) 


by  Lagrange's  formula  in  the  form 


au.  n  .  au. 

— i  An'*  +  —L  a"’1 

an  1  an  2  ’ 


where  the  derivatives  (j  -  1,  2)  are  taken  at  some  intermediate  points 


nni=nr1+eni(nr1-r2)'  (i=i>2>> 


one  has  that 


I 
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•-*  tj* 


i  :tui 


mr^ 


•r-n  tr  ' 


..2  n  n-2.2  .  n-1  2  * 

M  max  t(n.  )  ,  (n  )  t 


^n  <  t _ 0<  t<  1 


An-1  , 
A  + 


.,2  /  n-2.2  .,2  .  n-1.2 

_  M  max  ( n.  )  _  M  max  ( n  ) 

+ 1 _ P-4SU _ A"-‘  + 1 _ 2^1 _ An 

2  ,5  2  ,5  ’ 

min  r  min  r  . 

0<t  <  1  n~2  0<t  <  1  n_1 


From  this,  it  follows  that 


..2  {.  n-2  2  .  n-1  2,  ..2  ,  n-2.2 

M  max  i(n.  )  ,(n  )  }  M  max  (£  ) 

T  0<.t<  1 _ +  T _ 0<  t<.l _  ^n-1 

2  min  fr5  r5  ,}  2  min  r5 

n  0<  t<  1  n-2  n-‘ _ 0<t<l  n~2  „ 

A  -  . .2  ,  n-1  2 

M  max  ( n  ) 

.  T _ 0<.t<  1 _ 

1  "  2  ,5 

min  r  , 

0<  t<  1 


2.  .2  //  n-2,  2  /  n-1,  2t 

t  M  max  l(ji  )  ,  ( n  )  ) 

0<t  <1 _ 

min  (r5  r5  .} 

O^til  n-2  n-‘ 


2 

‘-T 


2  .  n-1  2 

M  max  (ri  ) 

0<  t<  1 


min  r 
0<  t  <  1  1 


An_1< 


2^2  /  M 

t  M  max  ( n.  ) 
0<t<l _ 

min 

0<t<l  k 

M  max  (n  ) 

T  0<Jt^l 

2  2  5 
min  r, 

0<  t<  1  k 


A  (24) 


k  =  0, 1,  2, . . . ,  n  -  1  , 


under  the  hypothesis  that 


.  .2  k  2 

M  max  ( ri  ) 

0<  t  <  1 


<  1  . 


0  <  t<  1 
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*»  x  >, 

'  ' a  . rv 


We  shall  prove  by  mathematical  induction  that  every 
rn  _>  kh  and  nn  <_  k1  n_°  , 

where  k  is  a  positive  number  less  than  1,  and  k1  is  a  positive  number  greater 


than  1 . 


For  n  =  0,  these  inequalities  are  obviously  valid,  in  so  far  as  r^  >  h. 


We  shall  assume  that  it  is  valid  for  all  r^  and  n^  for  j  =  1,  2, . . . ,  n-1, 


so  that 


rj  ^  kh  >  fVH0  ( j  =  1,  2, . . . ,  n-1)  . 


(26) 


If  it  follows  from  this  that  the  inequality  (  26)  is  satisfied  for  j  =  n,  then 
it  will  be  proved  for  all  n. 


By  hypothesis, 


2  »*2,  1 2  0 
t  M  k1  n 


2  A5 

k  h 


by  virtue  of  (  24)  and  (  26) 


2./I2,  .2  02 

t  M  k1  nu 


A5 

k  h 


.  T2  lA-V  I 
'  2  kV  y 


where,  as  can  be  seen  from  equation  ( 16), 

,  2  0 

1  t  ii 


1  1  tZM2u! 
1  ~~2  5 

ro 
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o 


»/ 


and  rQ  and  n_  satisfy  the  condition 


T2M2n° 

2r! 


<1  . 


Equation  (  23)  and  relationships  (  26) ,  (  28)  ,  (  29)  allow  one  to  write  the 

following  inequality: 

2  n-1  2,  2  .  n-1.  2 

t  max  n_  t  M  max  (n  ) 

lnn-n°l  li  - +  4 - l- 

min  r  .  min  r  k  h 

0<  t  <  1  n"  0<  t<  1  n_1 


f  T2M2k'2H0 


n-1 


1  22,  ,2  0 

_1_  t  M  k1  n 

2  iA5 
k  h 


.5,5 
k  h 


V 


2  * .  2.  ,2n0 

t  M  k1  H 

1  "  2  5  5 


2 


k  h  / 


.  2,  ,  0 

_1  t  k1  n_ 

8  A3 

k  h 


1  - 


i2M2n° 


It  is  required  that 


2  2  . 2  0 
t  M  k1  H 

,5.5 
k  h 


1- 


t2  M2k'2n° 


,5.5 
k  h 


<  1  , 


or 


,  2  2  . 2  0 
3_  t  M  k*  XL 

2  A5 

k  h 


<  1  . 


(30) 


■s  . 


Note  that  inequality  (  30)  guarantees  the  satisfaction  of 


I  n  Oil  -r^k'n0 

in  -n  I <s  -rr 

k  h 


1  +  i 


1  - 


1  T2M2n° 


27).  Now, 
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v~ 


and 


n  .  0  ,  i  T^k'J!0 
n  <n  +- 


8  3  3 

k  h 


1+T 


l 


2  2  0 
1  T  M  H 

1  “  2  5 


0 


We  obtain  an  estimate  for  r  , 

n  * 


n  .  n 


.....  .  n  0.  ,  ,  n  0,  ,0  ,0 

•En=Vl  +  n2~r2  s«ni-ni)£i  +  <n2-n2>-rz  +  niil+n2X2“ 


,  .  n  o,  .  n  0k 

=J0+,n1  -n1)X,  +  (n2-n2)r2  , 


Ir  -r  I  <1  nn-  n°  I 
~  n  ~  0 - — 


r,+  rj 


1  r2k,n° 


1  2  —  8 


3  3 
k  h 


1+7 


1  T2M2n° 


(rl  +  r2> 


and 


.  2.,  0 

^  1  t  k'a 

n-  0  8  33 

k  h 


1+T 


1  T2M2n° 


1  - 


<r,  +  r2)  . 


n 


In  so  far  as  rQ  >  h , 


r  >  h  - 


1  T2k»n° 


n-  8 


3  3 
k  h 


1  +T 


1 


1  t2M2U° 
1  "  2  5 


(  rl  +  r2> 


0 


* 


: 


We  determine  k  and  k'  so  that 


2  0 
0  ^  1  t  k'nu 

—  +8  ,  3,  3 


kh' 


*4 


1  -4 


1  tV/ 


<  k*  n 


0 
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5 


t. 


* 


I 


n  0 

Then  in  ^  k1 11  ,  r^^kh,  which  proves  inequality  (  26)  for  arbitrary  n. 
Inequalities  (  31)  and  (  32)  may  be  written  in  the  following  way: 


k'  > 


1 


8  .3.3 
k  h 


1+T 


l 


,  2»,2  0 

.  _1_  t  M  XL 

1  "  2  5 


‘0 


.  2,  ,0 

k  <  1  -±LJ£L- 

-  4  .3.  3 

k  h 


>♦* 


1  - 


1  x2M2n° 


(  33) 


(34) 


l 

1 

1 


) 


If  k  and  k'  may  be  chosen  in  order  that  (  30) ,  (  33),  (  34)  are  satisfied, 
then  the  process  of  successive  approximations  defined  by  formulas  (  20)  converges. 

From  (  33)  and  (  34) ,  assuming  that  equality  holds,  k  and  k'  may  be 
determined  by  the  method  of  successive  approximations,  taking  k  =  k*  =  1  as 
initial  approximations. 

It  is  possible  to  obtain  convergence  conditions  more  simply,  which,  however, 
are  cruder. 

If  one  chooses  initially  that 


1  t  Vn° 


) 


then  for  this,  the  satisfaction  of  inequalities  (  33),  (  34)  is  sufficient  to 
guarantee  the  satisfaction  of  the  following  inequalities: 
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T 


! 


k'  > 


1  - 


5  T 

24  ,3,3 
k  h 


k  < 


1  - 


5  t  2  k 1  n° 
12  k3h3 


Setting  k  =  j  ,  k'  =  —  ,  the  convergence  conditions  take  the  form 


h3  >  10t2 


h°  >  108M2t2  . 


I 


By  virtue  of  inequality  (  24) ,  if  condition  (  30)  is  satisfied,  the  sequences 
of  functions  In^1},  (n^l  converge  to  functions  n^n^,  which  are  solutions  of 
the  system  (  18) . 

One  may  convince  oneself  of  this  immediately  by  letting  n  - 00  in 
equations  (  20) . 

The  legality  of  passing  to  the  limit  under  the  integral  sign  follows  from  the 
uniform  convergence  of  the  functions  in  the  integrand,  in  as  much  as  the  estimate 
(  24)  is  independent  of  time. 

We  shall  prove  the  uniqueness  of  the  solution  obtained  in  the  region 


3_ 

2 


2  2  2 
t  M  max  n. 

0£t<l 

min  r^ 

0<t  <  1 


<  1  . 


Suppose  that  there  are  other  solutions  n| ,  n’^  which  satisfy  condition  (  30) 
Consider  the  differences 
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In  as  much  as  the  expression  in  parentheses  is  positive,  one  obtains 

I  n_'  -  nj  <_  0  , 

which,  by  the  virtue  of  the  non-negativity  of  |n_'  -  nj  is  only  possible  if 

I  n*  -  n|  =  0  . 
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In  order  to  obtain  an  estimate  for  the  error  committed  if  the  calculation  is 


stopped  at  the  nth  approximation,  we  consider  the  difference 


00  .  n-1  . 

k  a  .1 


I  n  -  nn|  =max  (1  nj- n^| ,  |  n2-n2!  }_<  2  A  =  A  , 

k  =n 


where 


a  = 


2  2  l2n° 

t  M  k1  H 

.5.5 
k  h 


2  2  .2  0 

t_  M  k1  H 

1  "  2  .  5.  5 


k  h 


x  and  k'  satisfy  inequalities  (  33) ,  (  34),  h  satisfies  (  30) ,  and  A'  satisfies 
inequality  (  29) . 


§5.  The  Case  of  Perturbed  Motion. 

If  in  addition  to  the  attraction  of  the  central  body,  a  force  mF  acts  on 
the  planet  or  satellite,  then,  instead  of  equation  (1),  the  following  equation 
of  motion  holds: 

2 

r(t)  =  k  ( 1+  m)— r -  +  F(  t)  .  (35) 

r  ( t) 


The  corresponding  integral  equation  may  be  written  in  the  form 


21  rr(t')  . 

r(t)=r  (t)  +  x  Jk(  t,  t')(  — -  -F(fidf  . 

0  W(t')  1 


(36) 


If  the  force  mF  is  the  perturbing  action  of  a  planet  with  mass  m  and 

s 

radius  vector  r  ,  equation  (  36)  takes  the  form 
**  s 
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2  1  j  r(t')  m 


r(f)-r  (f)  r  (tf) 
s  ■+■  s 


Llr(f)-r  (f)|J  rV)J 


> dt1  .  (  37) 


Set 


r(  t)  =  n  ( t)  r  +  n  ( t)r  +n  ( t)r 
~  1  ~  1  2  ~  2  s  ~s 


then  equation  (  37)  reduces  to  a  system  of  three  equations  for  the  quantities 


n!<  t) ,  n2(  t) ,  ng{  t)  : 


n}(  t)  =  n°(  t)  +  t2 


1  n(f) 

I  K(t,t')-^ -  1  + 

6  r  (f) L 


m  3. ... 

s  r  ( t1) 


n2(t)  =n2(t)  +  r2 


>  *  n2(t)r 

0  r  ( t')  L 


(f)  L  1+m  I  r  ( t* )  -  r  (t')l3 
~  ~  s 

n’(  ^  r  r3(f) 


1  +■ 


(  f)  L  1+m  |r(t')-r  (t')|3  J 

~  ~  s 


dt' 


dt' 


1  fn(t')  mc  n  ( t' )  —  1 


n  (t)  =  r  J  K(  t,  t'  K  ”7 - + 


2r‘ . I  ‘V^  .  ‘“s  “s'1''-1  .  ms  1 


r3(f)  1+m 


I  r  ( t' )  -  r  (t')|  1+m  r  (t')J 

~  ^  s  s 


)df  , 


which  may  be  solved  by  Newton's  method. 


0  0 


As  initial  approximations  one  may  take  n.  ,n,,n  =0  as  in  the  case  of 

1  2 '  s 

unperturbed  motion. 


§6.  Determination  of  Orbits  from  Three  Observations 
Let  tj,  t^,  t^  ( tj  ^  t2  _<  t^)  denote  the  times  of  three  observations  chosen 
for  the  determination  of  an  orbit.  Each  of  the  three  observations  permits  writing 
a  relationship 

~i  "pi~i  -Sj  (‘*1,2,3),  (38) 

where  p  e  ,  denotes  the  geocentric  radius  vector  of  the  planet  (  here,  e  is 

i  ~  i  7  ~  i 

the  unit  vector  indicating  the  direction  to  the  planet,  pi  is  the  geocentric  distance) , 
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and  R  is  the  geocentric  solar  radius  vector. 

~  i 

For  the  determination  of  the  orbit  given  the  directions  e  to  the  planet 
and  the  geocentric  solar  radius  vectors  R  ( i  =1,2,3)  the  boundary  values 
r  and  r,  may  themselves  be  unknowns,  and  have  to  be  determined  from  some 
system  of  equations. 

To  obtain  these  equations,  one  uses  the  relationships  which  connect  the 
geocentric  coordinates  for  three  given  times 


r  =  n  ( t  )  r  +  n  ( t  )  r 
~  2  1  2  ~  1  3  2  ~  3 


or,  by  virtue  of  (  38) , 


p2~  2  ”  ~  2  "  pl-  l“  5  +  n3^  V  ^  p2~  2~  ~  2^  * 

[  e  e  I  [e.e  .] 

Multiplying  this  equation  first  by  ---  — — - — —  ,  then  by 


e  .[  e  e  ] 
-1  ~ 


e.fe  e  ]  ’ 
~  1  ~  2  ~  3 


-  ,  one 


obtains 


nl(t2)  P1  "  ^  ni(  V5l  +  n3^  t2)~3_  5  2^^6263]  » 

n3(  t2)p2  ni(  t2)  ^1  +  n3(  t2)  ~  3  ~  ~  2^  e.[e  e  ] 

~  1  ~  2  ~3J 


(39) 


e  .[  e  e  ]  *  0  . 
~  r  ~  2~  3 


The  solution  of  system  (39),  together  with  the  system 


Uj(  t)  =  n°(  t)  +  T^  /K(  t,  t1) 


nj(  f) 


nj(  t'>i1+n2(  t*)x2l 


dt*  , 


n3(  t)  =n3(  t)  +r  f  K(  t,  t') 


n3(f) 


0  1  t')^1+n3(  tf)^3 


“dt'  O 


(40) 
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determines  and  p^  ,  and,  consequently,  and  r^. 

As  the  initial  approximation  for  the  ratio  of  the  area  of  the  triangles,  one 
chooses  n^,n^,  which  are  ratios  of  intermediate  times 


0 


n 


1 


where 

=  Wl  +  m  ( t3  -  t2)  , 
t2  *  Wl  +  m  (t3  -  tj)  , 
T3=kv/l  +  m(t2-t  )  . 


Substituting  them  into  equation  (  39),  one  obtains  p^,  p^  in  the  first  approxi¬ 


(  38) ,  and  equations  (  40)  are  solved  with  them  to  obtain  new  values  of  n^,  n^, 
and  so  on.  Therefore,  equations  (  38)  -(  40)  together  permit  the  determination 
of  the  geocentric  distances  p^,  p^,  p^,  and,  consequently,  the  geocentric  radius 
vectors  £^>£2,~3  t^ie  P^anet-  Knowing  these  values,  one  may  determine 
the  elements  of  the  orbit  in  the  ordinary  way. 


mation. 

Approximate  values  of  are  determined  with  the  aid  of  relationships 
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